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DIFFUSION  THEORY  FOR  ADSORPTION  AND  DESORPTION  OF  GAS  ATOMS  AT  SURFACES 

* 

Patrick  J.  Pagnl  and  James  C.  Keck 
Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts 


ABSTRACT 

'The  time  evolution  of  a  nonequilibrium  ensemble  of  gas  atoms  adsorbed 
on  a  solid  surface  is  described.  Adsorption  and  desorption  rate  constants 
and  adatom  (adsorbed  atom)  energy  distributions  are  obtained  in  the  steady 
state  approximation  from  the  appropriate  master  equation  or  an  equivalent 
diffusion  equation.  Energy  transition  probabilities  are  obtained  from  a 
classical  gas-surface  collision  model  using  an  interaction  potential  com¬ 
posed  of  an  oscillating  harmonic  repulsion  and  a  stationary  attraction.' 

Calculations  have  been  performed  over  the  following  range  of  similarity 

_2 

parameters:  inertia  ratio  10  <  p  =  m  /m  <2,  frequency  ratio 

g  s 

1  <  v  ■  to  /to  <  10,  and  well  depth  2  <  6  ■  D/kT  <  100.  For  "resonance" 
s  g 

values  of  the  frequency  ratio,  gas-surface  collisions  are  sufficiently 
adiabatic  that  adatom  energy  distributions  are  significantly  nonequilibrium, 
and  steady  state  rate  constants  are  depressed  below  the  equilibrium  rate 
constants.  Applications  include  prediction  of  adsorption  and  desorption 
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rate  constants,  sticking  and  accommodation  coefficients,  sublimation 
and  condensation  rates,  and  boundary  conditions  for  the  Knudsen  layer. 

I .  INTRODUCTION 

This  investigation  concerns  the  evolution  of  a  system  of  gas  atoms 
adsorbed  on  a  solid  surface  when  that  system  is  thrown  out  of  equilibrium. 

A  fundamental  physical  understanding  of  the  nonequilibrium  phenomena  at 
such  a  surface  is  needed  to  delineate  important  similarity  parameters  and 
provide  guidelines  for  definitive  experiments  in  surface  studies.  Consider, 
as  an  initial  problem  statement,  an  ensemble  of  gas  atoms  in  equilibrium 
with  a  solid  surface  at  temperature  T  enclosed  in  a  sealed  container.  The 
enclosure  is  instantaneously  evacuated  at  t  »  0.  The  problem  is  to  deter¬ 
mine  the  time  history  of  the  surface  concentration  and  the  adatom  (adsorbed 
atom)  energy  distribution.  Experiments  show  that  over  a  wide  range  of 
conditions  the  desorption  process  proceeds  according  to  a  first-order  rate 
equation: 

k  F(t)  -  k.  A(t)  (1) 

dt  a  d 

_2 

where  A  is  the  surface  concentration  (cm  ),  F  is  the  gas  concentration 
(cm  ),  and  kg  and  k^  are  the  adsorption  and  desorption  rate  constants. 

A  satisfactory  quantitative  understanding  of  temperature  and  atomic 
species  dependence  of  the  rate  constants  kg  and  is  sought.  This 
analysis  is  the  first  attempt  the  authors  are  aware  of  to  describe  the 
functional  dependence  of  and  k^  using  nonequilihrium  kinetics.  An 
equilibrium  approach^  which  assuntps  a  Roitzrann  dfstribution  for  the 
adatom  energy  distribution  and  that  all  atoms  which  strike  the  surface 
adsorb  gives  the  classical  desorption  rate  constant 
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where  w  is  a  characteristic  frequency,  D  is  the  potential  energy  well 
8 

depth,  and  T  is  the  surface  temperature.  The  statistical  description 
given  here  approaches  this  equilibrium  expression  as  a  limiting 
case. 

Many  previous  investigations  of  the  interactions  of  gas  atoms  at 

(2-21) 

surfaces  have  been  reported,  and  excellent  reviews  of  this  work 

(22)  (23) 

may  be  found  in  articles  by  Trilling  and  Kogan.  Of  particular 

interest  in  connection  with  the  present  investigation  is  the  work  of 
(2-9) 

Armand  who  has  made  careful  studies  of  the  effect  of  lattice  interac- 

(18) 

tions  on  the  desorption  process  and  Logan  and  Keckv  1  who  have  developed 
a  relatively  simple  "soft  cube"  model  for  scattering  of  gas  atoms  by 
surfaces.  Also  of  interest  is  the  analogous  gas-phase  problem  of  molecu- 

/  ry  f  O  1  \ 

lar  dissociation  and  recombination,  and  a  review  of  this  litera- 

(32) 

ture  has  been  published  by  Bunker.  The  theory  developed  in  this 

(27) 

paper  is  based  on  this  analogy  and  parallels  the  work  of  Keck  and  Carrier 
and  Keck.(28) 

Applications  of  the  theory  include  prediction  and  correlation  of 
adsorption  and  desorption  lifetimes,  accommodation  and  sticking  coeffi¬ 
cients,  sublimation  and  condensation  rates,  and  scattering  patterns.  In 

addition,  potential  parameters  for  gas-surface  interactions  may  be  deduced 

(33-39) 

from  a  comparison  of  the  theoretical  results  with  experimental  data 
obtained  using  modulated  beam  and  flash  desorption  techniques. 

A  general  discussion  of  the  theoretical  model  and  the  governing 
master  equation  used  to  obtain  rate  constants  is  given  in  the  following 
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secclon.  The  mechanism  of  the  gas-surface  collisions  aid  the  dimension¬ 
less  parameters  characterizing  the  problem  are  discussed  in  Section  III. 
The  transition  kernel  appearing  in  the  master  equation  and  the  steady 
state  adsorption  and  desorption  rate  constants  are  derived  in  Section  IV. 
Finally,  a  brief  summary  and  conclusions  are  given  in  the  last  section. 

II.  MATHEMATICAL  MODEL 
A.  Qualitative  Discussion 

Consider  first  a  single  gas  atom  approaching  the  solid  surface  from 
the  free  state  with  a  nondimensional  energy  e  -  E/kT  as  shown  in  Figure  1 
where  r  is  the  atom  displacement  from  the  surface.  Far  from  the  surface, 
the  free  atom  feels  an  attractive  van  der  Waals  force  represented  as  a 
potential  with  a  positive  gradient  increasing  as  the  separation  distance 
decreases.  In  the  Immediate  vicinity  of  the  surface,  the  gas  atom 
encounters  a  repulsive  valence  force  represented  by  a  potential  with  a 
negative  gradient  that  approaches  infinity  as  the  separation  distance 
goes  to  zero.  At  finite  separation  distances  these  two  poentlals  over¬ 
lap  to  form  a  potential  energy  well  which  defines  the  collision  mechanics. 
For  the  collision  shown,  the  atom's  energy  after  collision,  e',  is  less 
than  the  well  depth,  and  the  atom  is  trapped  on  the  surface. 

The  distinctive  feature  of  the  collision  between  a  free  atom  and  a 
lattice  atom  is  that  the  "inner"  potential  oscillates  while  the  "outer" 
potential  is  stationary.  The  stationary  outer  potential  results  from 
averaging  over  pairwise  potentials  between  the  gas  atom  and  several 
lattice  atoms  which,  at  large  separation  distances,  are  equidistant  from 
the  gas  atom  and  hence  have  an  equal  influence  on  it.  The  oscillations 
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of  the  repulsive  potential  are  not  eliminated  by  averaging,  since  at 
small  separation  distances  one  lattice  atom  dominates  the  repulsion. 

Now  consider  an  ensemble  of  adatoms  trapped  on  the  surface.  The 
array  of  surface  atoms  is  modeled  as  a  set  of  independent  simple  harmonic 
oscillators  and  may  be  visualized  as  a  "simple  harmonic  chess  board"  in 
which  each  square  oscillates  perpendicular  to  the  plane  of  the  board 
with  a  frequency  approximated  by  a  surface  Debye  temperature  and  with 
an  amplitude  probability  specified  by  a  Boltzmann  energy  distribution 
at  the  temperature  of  the  solid.  The  well  depth,  which  is  a  function 
of  both  time  and  position  on  the  surface,  is  approximated  by  a  constant. 
This  approximation  is  possible  since  all  important  energy  levels  are 
near  the  top  of  the  well.  A  free  atom  at  rest  an  infinite  distance 
from  the  surface  is  assigned  an  energy  equal  to  the  nondimensional  well 
depth  6  ■  D/kT,  placing  the  zero  of  energy  at  the  bottom  of  the  well. 

The  adatom  energy  distribution  can  be  found  from  a  kinetic  master  equa¬ 
tion  once  energy  transition  probabilities  are  obtained  from  a  gao-suv- 
face  collision  model.  This  is  equivalent  to  treating  the  movement  of 
adatoms  among  the  energy  levels  of  the  potential  well  as  a  random  walk 
process.  Adatom  energy  transitions  are  statistically  independent  since 

the  relaxation  time  of  a  struck  lattice  atom,  T  ,  and  the  collision 

>  r> 

duration,  Tc>  are  less  than  the  period,  Tq,  of  the  adatom  oscillation 
in  the  gas-surface  potential  well.  The  kinetic  equation  used  here  only 
accounts  for  adatom-surface  collision  Induced  energy  transitions,  thus 
this  analysis  is  restricted  to  low  surface  concentrations. 

The  desorption  relaxation  process  has  three  distinct  time  scales. 
Initially  there  is  a  transient  regime  of  the  order  of  the  characteristic 
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adatom  oscillation  time,  T  ,  during  which  there  is  negligible  desorption. 

During  the  transient,  the  adatom  energy  distribution  relaxes  to  a  steady 

state  form  which  is  the  asymptotic  energy  distribution  as  t/TQ  •*  °°.  The 

second  time  regime  is  characterized  by  the  mean  desorption  lifetime, 

T,  =  k  *.  Since  T,  is  of  the  order  of  T  exp(D/kT)  and  the  well  depth 
d  d  d  o 

is  large  compared  to  kT,  »  T  .  This  difference  in  time  scales  allows 
the  steady  state  distribution  to  be  used  during  the  entire  second  time 
regime  while  all  significant  desorption  takes  place.  The  third  regime  is 
the  equilibrium  condition  at  t/T^  -*■  °°.  Here  the  adsorbing  flux  equals  the 
desorbing  flux.  It  is  the  second  time  regime  which  will  be  examined  in 
detail  to  determine  the  temperature  and  atomic  species  dependence  of  the 
rate  constants. 

B.  Master  Equation 

Assuming  that  the  state  of  an  adatom  is  specified  by  its  energy, 
the  kinetic  master  equation  for  the  reduced  adatom  energy  population 
X(e,  t)  can  be  written: 

Ne(e)  fj-  X(e,  t)  -  /q  R(e,  e’)  [X(e\  t)  -  X(e,  t)]de' 

+  R(e,  f)  [H(e,  t)  -  X(e,  t)]  (3) 

where 

X(e,  t)  -  N(e,  t)/Ne(e)  (4) 

is  the  ratio  of  the  actual  energy  distribution,  N(e,  t),  of  adsorbed  atoms 
to  an  equilibrium  distribution,  Ne(e)  =  (N/Q)e  e(jifl/de),  at  the  surface 
temperature  T,  N  is  the  number  of  atoms  in  that  region  of  phase  space  for 
which  the  partition  is  Q,  dStyde  is  the  energy  level  density  per  unit  e. 
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e  =  E/kT  is  the  nondimens ional  adatom  energy,  6  ■  D/kT  is  the  nondimen- 
sional  well  depth,  R(e,  e'Jdede*  is  the  one-way  equilibrium  transition 
rate  from  an  energy  band  de  at  e  to  an  energy  band  de'  at  e', 

R(e,  f)de  *  R(e,  e')de'de  (5) 

is  the  one-way  equilibrium  desorption  rate  from  an  energy  band  de  at 
e  to  the  free  state,  and  H(e,  t)  is  the  ratio  of  the  actual  one-way 
adsorption  rate  to  the  equilibrium  one-way  adsorption  rate,  R(f,  e). 

Note  that  by  detailed  balancing,  we  have  R(e,  e')  =  R(e',  e)  and 
R(f,  e)  -  R(e,  f). 

X(e,  t)  is  a  useful  measure  of  nonequilibrium  since  it  is  unity 
at  equilibrium.  X(e,  t)  >  1  indicates  over  population  of  the  energy 
level  e,  and  X(e,  t)  <  1  indicates  depletion  of  the  level  e.  The 
master  equation  simply  states  that  the  adatoms  leave  the  energy  level 
e  either  by  making  transitions  to  the  level  e*  at  a  rate  R(e,  e')  X(e,  t) 
or  by  desorbing  to  the  free  state  at  a  rate  R(e,  f)  X(e,  t)  and  enter 
the  energy  level  e  either  by  making  transitions  from  the  level  e'  at  a 
rate  R(e',  e)  X(e',  t)  or  by  adsorbing  from  the  free  state  at  a  rate 
R(f,  e)  H(e,  t). 

The  form  of  Equationf 3) results  from  the  following  assumptions: 

1.  The  adatom  concentration  is  sufficiently  small  that  collisions 

between  adatoms  may  be  neglected.  If  an  equivalent  Knudsen  number 

Kn,  is  defined  as  the  ratio  of  the  adatom-surface  collision  rate 

to  the  adatom-adatom  collision  rate,  this  assumption  is  valid  when 

Kn  »  1.  The  adatom-surface  collision  rate  is  the  inverse  period 

13  -1 

of  oscillation  in  the  well '>'10  sec  .  The  adatom-adatom  colli¬ 
sion  rate  is  a  c  A  where  A  is  the  surface  concentration,  a  the 
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collision  diameter  <\<  10  cm,  and  the  mean  speed  c  ■  (irkT/2m  'ulO"* 

8 

cm  sec.  The  effective  Knudsen  number  is  then  Kn  -  T  ~^/(o  c  A)  -  10^  fTl 

o 

16  2 

and  assumption  1  is  valid  if  A  «  10  cm”  . 

2.  A  classical  approach  is  valid.  Keck  and  Carrier have  shown  that 
quantum  models  fail  when  the  average  energy  transfer  per  collision 
is  greater  than  the  energy  level  spacing  at  the  rate-limiting  step 
in  the  well.  For  the  transition  rate  derived  here,  the  "bottleneck" 
is  at  the  top  of  the  well  where  the  energy  level  density  is  high, 

and  the  quantum  number  is  large.  For  very  shallow  well  depths,  6  %  1, 
this  approximation  is  not  valid. 

3.  The  nonequilibrium  transition  probability  is  equal  to  the  equilibrium 
transition  probability.  This  assumption  is  always  valid  since  the 
collision  mechanics  are  determined  by  the  temperature  of  the  solid 
and  the  atomic  parameters  of  the  gas  and  solid  and  are  independent 

of  the  energy  distribution  of  the  adsorbed  ensemble. 

C.  Steady  State  Solutions 

To  solve  Equation  (3)we  now  Introduce  the  steady  state  approximation 
by  setting  the  time  derivative  on  the  right-hand  side  equal  to  zero. 

This  gives 

/£  R(e,  e*)  [X(e * ,  t)  -  X(e,  t) Jdc*  +  R(e,  f)  [H(e,  t)  -  X(e,  t)J  -  0  .  (7) 

The  form  of  Equation  (7j  Indicates  that  energy  and  tiu'e  variables  will 
separate  if  it  is  assumed  that  a  steady  state  reduced  population,  x(e)» 
exists  such  that 

x(e)  f(t)  -  X(c,  t)  -  H(e,  t)  :  e  <  6  (8) 

where  f(t)  gives  the  time  dependence  of  the  adsorbed  population.  If 
the  energy  distribution  in  the  free  gas  is  steady,  H(e,  t)  may  also  be 
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separated  as  H(e)  F (t) /Fe  where  F(t)  is  the  time-dependent  free  gas 
concentration,  and  H(e)  is  the  ratio  of  the  adsorption  rate  at  an 
energy  level  e  from  the  actual  incident  distribution  to  the  rate  from 
an  equilibrium  distribution.  In  terms  of  the  equilibrium  transition 
rate  R(£',  e)  and  a  time-independent  incident  atom  energy  distribu¬ 
tion  xg(e*)  =  Fe  Ng(e',t)/F(t)  Ne(e’), 

H(e)  -  /“ R(e',  e)  X  (e')  de'/R(f,  e)  .  (9) 

H(e)  is  unity  when  the  incident  energy  distribution  is  Boltzman  about 

the  temperature  of  the  solid  so  that  X(e')  *  1. 

g 

An  integral  equation  for  X(e)  is  obtained  by  substituting  Equation 
(8) into  Equation (7) and  assuming  9H/3e  «  f(t)  dX/de 

R(e,  e')  (x(e')  -  X(e))  de’  +  R(e,  f)  x(e)  -  0  .  (10) 

The  condition  3ll/3e  «  f(t)  dX/de  is  generally  satisfied  since  near  the 
top  of  the  well  the  gradient  of  X(e)  is  large  due  to  depletion  of  the 
upper  energy  levels  by  desorbing  adatoms, and  11(e)  does  not  have  a  large 
gradient  at  any  energy.  The  steady  state  population,  X(e),  is  obtained 
from  a  iteration  solution  to  EquationflO^ or  by  conversion  of  Equation (10) 
to  an  equivalent  diffusion  equation . 

The  surface  concentration  is  given  by 

A(t)  =  f{)  X(e,  t)  N0(e)  de  ,  (11) 

where  X(e,  t)  is  obtained  in  the  steady  state  approximation  by  assuming 
X(e,  t)  is  separable,  and  setting  X(0)  =  1  in  Equation^.  This  gives 

X(e,  t)  -  X(0,  t)  X(e)  +  (H(e)  -  H(0)  X(e))  E(t)/Fe  .  (12) 
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Sin  ce  Nfi(e)  Is  a  rapidly  decreasing  function  of  C,  substituting  Equa¬ 
tion  (12)  into  Equation  (11)  with  H(e)  slowly  varying  and  x(0  ~  1 
near  e  ■  0  gives 

A(t) ~  X(0,  t)  /q  Ne(e)  de  -  X(0,  t)Afi  .  (13) 

The  rate  constants  are  obtained  as  functions  of  transition  rates  R(f,  e), 
and  the  steady  state  reduced  population,  x(e)»  from  the  energy  Integral 
of  the  master  equation.  Integrating  the  original  master  equation  (3) 
over  energy  using  Equations  (9)  and  (13)  yields  the  phenomenological 
equation 


dA(t) 

dt 


-  kfl  F(t)  -  kd  A(t)  , 


where  the  rate  constants  have  the  form 


(14) 


ka  -  f60  R(f,  e)  Xfe)  d£/Fe  (15) 

and 

kd  -  /q  R(f ,  e)  xte)  de/Ae  06) 

if  H(e)  is  unity.  These  are  the  steady  state  rate  constants.  The 
equilibrium  rate  constants  are  given  by  setting  *  1  in  Equations 

(15)  and  (16).  It  is  seen  from  these  equations  that  the  ratio  of  the 
forward  to  backward  rates  in  both  the  steady  state  and  equilibrium  case 
is  equal  to  the  equilibrium  constant. 


F 

K  -  r - 
e  A 

e 


de  _  d 

k  k 
ae  a 


(17) 


where  Kfi  is  the  equilibrium  constant,  kdg  is  the  equilibrium  desorption 

rate  constant,  and  k  is  the  equilibrium  adsorption  rate  constant.  Whc 

ae 

H(e)  is  not  unity;  i.e.,  an  incident  gas  stream  exists  which  does  not 
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have  a  Boltzman  distribution  about  the  temperature  of  the  solid,  the 
ratio  of  the  forward  to  backward  rate  constants  is  equal  to  the  product 
of  the  equilibrium  constant  and  11(0), 

kd/ka  “  Kc  ,!<0)  *  (18) 

The  factor  11(0)  takes  into  account  the  dependence  of  the  adsorption  rate 

constant  on  the  energy  distribution  of  the  incident  atoms.  In  this 
case,  the  energy  integral  of  the  original  master  equation  with  Equa¬ 
tions  (12),  (13),  and  (14)  gives  the  adsorption  rate  constant  as 

ka  "^0  R(f*  G)  X(C>  dG)  H<°>/Fc  •  (19) 

Since  the  desorption  rate  constant  is  unaffected  by  the  incident  atom 
energy  distribution,  it  remains  unchanged  and  is  given  by  Equation  (16), 

III.  COLLISION  MECHANICS 
A.  Equations  of  Motion 

The  mechanics  of  a  single  collision  between  a  gas  atom  and  a  sur¬ 
face  atom  must  be  described  in  detail  to  obtain  the  energy  transition 

probability,  R(c',  c).  The  collision  model  used  is  similar  to  the  soft 

(18) 

cube  model  of  Logan  and  K«2k  with  an  oscillating  half  harmonic 
repulsion  replacing  the  oscillating  exponential  repulsion  since  lower 
incident  energies  are  of  interest.  It  is  assumed  that: 

1.  tangential  momentum  is  conserved;  i.e.,  the  solid  surface  is  flat. 
This  is  not  true  since  the  crystal  structure  produces  a  rough  sur¬ 
face  on  the  atomic  scale.  However,  the  roughness  seen  by  an  inci¬ 
dent  atom  is  energy  dependent and  at  incident  energies  <  5  the 
adatom,  remaining  relatively  far  outside  the  crystal  structure, 
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would  see  an  approximately  flat  surface.  Experimental  data  on 
out-of-plane  scattering  and  the  success  of  the  soft^)  and 
hard^^  cube  models  in  predicting  scattering  patterns 
support  this  assumption. 

2.  the  surface  atoms  are  independent  simple  harmonic  oscillators 

with  a  single  characteristic  frequency  obtained  by  approximating 

the  lattice  frequency  spectrum  with  a  delta  function  at  the  Debye 

frequency.  Decoupling  the  surface  array  into  independent  oscillators 

has  been  shown  to  have  small  effect  on  the  energy  exchange  process. 

It  is  also  assumed  that  only  a  single  surface  atom  is  involved  in  any 

given  collision.  It  has  been  shown  that  most  of  the  load  for  most  of 

(18) 

the  collisions  is  taken  by  a  single  surface  atom. v  '  When  the  adatom 
velocity  parallel  to  the  surface  is  so  large  that  the  distance  traveled 
during  a  collision  is  greater  than  the  lattice  space,  this  one-on- 
one  approximation  falls. 

3.  the  gas-surface  potential  consists  of  a  stationary  arbitrarily  shaped 
attraction  and  an  oscillating  quadratic  repulsion.  The  motivation 
for  this  assumption  is  discussed  in  the  introduction.  The  interac¬ 
tion  is  described  classically  and  would  not  strictly  apply  to  inci¬ 
dent  atoms  whose  DeBroglio  wavelength  is  of  the  order  of  the  lattice 
spacing  or  to  solid  temperatures  below  the  surface  Debye  temperature. 

The  model  interaction  potential  is  illustrated  in  Figure  2.  The 
inner  half  (r  <  0)  is  an  harmonic  repulsion  while  the  outer  half 
(r  >  0) is  a  static  barrier  of  height.  Y(t)  and  Z(t)  are  the  gas  and 
surface  atom  displacements  from  their  respective  equilibrium  positions 
which  are  separated  by  the  equilibrium  extension  of  the  gas-surface 
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"spring."  k  and  k  are  the  lattice  and  gas-ourface  spring  constants  , 
s  g 

mg  and  mg  are  the  surface  and  gas  atom  masses,  and  6  is  the  desorption 
energy  divided  by  kT.  The  time  origin  marks  the  initiation  of  a  colli¬ 
sion.  For  t  <  0,  the  gas  atom  is  force  free,  and  the  surface  oscillator 
is  in  simple  harmonic  motion. 


Z(t)  -  ~  sin  (u»  t  +  <f>) 

CO  8 

8 

and  (20) 

Y(t)  ■  -  ut  +  sin  $ 
s 

where  u  is  the  incident  gas  atom  velocity  in  the  well,  ug  is  the  angular 
frequency  of  the  surface  atom,  and  V  and  <p  are  the  velocity  amplitude 
and  phase  of  the  surface  atom.  The  initial  conditions  are  established 
as  the  gas  atom  impacts  at  the  end  of  the  gas-surface  spring.  At  t  *  0, 


(21) 


Y(0)  -  Z(0)  -  L-  gin  <(> 

8 

Y (0)  ■  -u  ;  Z(0)  ■  V  cos  <p  . 

The  equations  of  motion  during  collision  (t  >  0)  arc 
mg  Y(t)  -  -  kg(Y(t)  -  Z(t)) 

m  Z(t)  -  -  k  Z(t)  +  k  (Y (t)  -  Z(t) )  , 

9  a  5 

The  collision  ends  when  the  gas-surface  spring  returns  to  its  equilibrium 

extension  with  the  gas  atom  moving  away  from  the  surface.  For  t  ■  T  , 

c 


(22) 


Y(T  )  -  Z(T  ) 
c  c 


(23) 


with  Y(lc)  positive  and  greater  than  Z(Tc).  The  collision  duration,  Tc, 
is  the  smallest  non-zero  root  of  Equation  (23).  The  nondimensional 
energy  transfer,  £(e),  is  defined  as  the  Increase  in  the  gas  atom  kinetic 
energy  during  the  collision  divided  by  kT. 
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A(e)  -  e'  “  e  “  fkT  (Y(TC>2  “  Y(0)2)  * 


(24) 


B.  Similarity  Parameters 

Defining  w  as  a  characteristic  time  and  uou  *  as  a  characteris- 
8  8 

tic  length,  the  nondimens ional  equations  of  motion  corresponding  to 


Equations  (22)  are 


v  y  +  (y  -  z)  -  0 
v2  (z  +  z)  -  y(y  -  z)  -  0 


(25) 


where  p  ■  mg/®s  *8  the  inertia  parameter,  V  ■  is  the  frequency 

parameter,  z  *  Z/uto  *  Is  the  nondlmenslonal  surface  atom  displacement, 

y  -  Y/uui  Is  the  nondlmenslonal  gas  atom  displacement,  and  time  Is 
8 

nondlmenslonallzed  with  respect  to  w  The  nondlmenslonal  Initial 

8 

conditions  are 

y (0)  -  z(0)  -  V 


y (0)  -  -1  and  z(0)  -  V, 


(26) 


where 


y  .  V.cos.  t  .Vjslnjt  (27) 

C  U  8  U 

are  respectively  the  nondlmenslonal  velocity  and  displacement  of  the 
surface  atom  at  Impact.  One  set  of  Initial  conditions  applies  to  all 
Incident  energies  since  the  Initial  gas  atom  velocity  no  longer  appears 
explicitly  in  Equations  (25)  and  (26).  Although  the  surface  oscillator 
frequency  is  assumed  to  be  given  by  the  bulk  Debye  temperature,  ws  ■ 
any  surface  oscillator  frequency  may  be  used  to  obtain  w  from  the  dimen- 

o 

slonless  ratio  V  -  w  /u  .  The  original  nine  variables  V,  $,  w  ,  w  ,  m  , 

8  g  8  g  S 

m  ,  D,  T,  and  u  are  thus  reduced  *o  three  parameters  v,  p,  6,  and  two 
8 

Initial  conditions  V  and  V  . 

c  s 
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From  comparison  with  experiment,  the  parametric  ranges  of  Interest 

-2  2 
are:  10  <  u  <  2,  1  <  v  <  10,  and  2  <  <S  <  10  .  In  the  adiabatic 

limit,  v  -*■  <»  ,  the  energy  transfer  goes  to  zero  9ince  the  collision 
duration  is  large,  and  the  incident  particle  turns  around  slowly  losing 
no  energy.  In  the  impulsive  or  hard  cube  limit  v  0,  and  the  colli¬ 
sion  duration  goes  to  zero.  Here  the  adiabacity  factor  which  is  propor¬ 
tional  to  exp(-  a)  T  )  -►  1,  and  the  energy  transfer  i9  a  maximum, 
s  c 

C.  Normal  Mode  Solutions 

General  solutions  to  the  equations  of  motion  during  collision. 

Equations  (25),  may  be  written  as  a  linear  superposition  of  normal 
modes : 

J  U)  -  £  a.q.(t)  and  z(t)  -  £  B.q .  (t)  (28) 

j  J  J  j  3  3 

where  J  ■  1,  2  Indexes  the  two  normal  modes,  q^ ,  of  the  system.  The 
following  algebraic  equations  defining  the  normal  mode  frequencies  are 
obtained  by  substitution  of  possible  solutions,  y  -  a  eiU,t  and  z  -  8  elu)t, 
into  Equations  (25): 

(1  -  vV)  a  -  8  ■  0 

not  -  (v2  +  u  -  v2oj2)8  -  0  .  (29) 

The  resulting  normal  mode  frequencies,  and  nondimensional  with 
respect  to  wg,  are 

Wl,2  "  <<v2  +  V  *•  1)  -  C  <v2  +  u  +  l)2  -  4v2)1/2)1/2/Av  .  (30) 

The  amplitude  coefficients  in  Equation  (28)  are  related  by  8j  “  (1  -  w2v2)aj. 
Applying  the  initial  conditions  given  by  Equations  (26)  to  Equations  (25) 
yields  the  gas  atom  trajectory 
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2  2 

c  V  1  +  Vc  -  u)^ 

y(t)  “  -  : - 1  cos  u).t  -  - — sin  oi-t 

1-5  1  (1  -  O^v  1 


2  2 

V  1  +  V£  -  oj^v 

+  : — --  COS  0)„t  + - sr - 9  sin  0)9t  , 

1-5  2  (1  -  0*1  u2vz 


(31) 


and  the  surface  atom  trajectory: 


z(t)  -  - 


(i  -  u>2v2)cvo 

1  ~  S 


cos  w^t  -  — 


(1  -  u>2v2)<l  +  Vc  -  u^v2) 


(1  -  OojJv2 


-  sin  o)jt 
(32) 


(1  -  o>2v2)V 

+  - cos  ^t  + 


(1  -  oj2v2)(1  +  Vc  -  oj^v2) 
(1  -  C)ui2uj2v2 


sin 


u^t 


2,  2 


where  Z,  m  u^/u)^.  The  nondimensional  collision  duration,  Tc,  is  the 
smallest  non-zero  root  of  y(Tc)  "  z(Tc).  Using  Equations  (31)  and  (32), 
the  implicit  expression  for  Tc  is 


cosfa^  +  ^)  ■  - 


l“2l 

i^r 


C  cos(uj2tc  4>2) 


(33) 


where  a  and  $  are  the  amplitudes  and  phases  of  the  normal  modes.  In 
terms  of  the  initial  conditions  and  similarity  parameters,  these  are: 


V2  C2  (1  +  V  -  v2oj2)2 
8  ^  ^ _  c  2 


KT 


4(1  -  O2  r  4(1  -  C)2  0)Jv4 


2  2  2  2 
r  (1  +  V  -  VUT 
s  +  _ _  c  1 


4(1  -  C)2  '  4(1  -  ?)2  u>*w2v4 


(34) 


and 


.  -(1  +  V  -  v2oj2) 

<t».  «  tan"  { - = - 5 — } 

V  uS.  v 
8/1 


A  _1  “(1  +  V  -  V2W?) 

*2  ■  £)  . 

V,“l  “2  v 


(35) 
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Once  the  collision  duration  is  known,  the  energy  transfer  in  a  collision 
is  obtained  directly  from 


A  -  e[y(Tc)  -  1] 


where  e  ■  m  u  / 2kT. 
g 


(36) 


D.  Collision  Duration  and  Energy  Transfer 
Contours  of  constant  collision  duration  and  energy  transfer  in 
initial  condition  space  are  given  in  Figures  3a  and  3b.  The  time  con¬ 
tours  in  Figure  3a  are  given  for  relatively  impulsive  parameters, 
y  -  V  *  1,  while  the  energy  contours  in  Figure  3b  are  for  V  ■  3.16  and 
y  «  .5.  There  is  a  "plateau"  in  the  V  ,  V  plane,  separated  from  the 
remainder  of  the  space  by  a  discontinuity  along  the  heavy  dashed  lines 

•  a 

in  Figures  3a  and  3b  given  implicitly  by  y(T  )  -  z(T  )  and  y(T  )  -  z(t  ). 

c  c  c  c 

This  discontinuity  follows  from  the  definition  of  as  the  smallest  non¬ 
zero  root  of  a  periodic  function.  The  effect  of  higher  v  is  to  Increase 
the  size  of  the  plateau  and  the  range  of  initial  conditions  for  which 
adiabatic  collisions  with  longer  duration  occur.  The  approximately 
radial,  straight-line  behavior  of  contours  of  constant  collision  time 

in  Figure  3a  follows  from  the  definitions  of  V  and  V  .  Radial  lines 

s  c 

in  V  ,  V  space  are  roughly  lines  of  constant  phase  angle,  <J>.  Figure  3a 

8  C 

illustrates  that  the  duration  of  the  collision  is  largely  dependent  on 
the  phase  angle  of  the  surface  atoms  at  impact.  Sln/ce  the  initial  condi¬ 
tions  establish  the  nature  of  the  collision,  regions  in  V  ,  V  space  may 

c  s 

be  assigned  a  particular  type  of  collision.  For  V  <  -  1,  region  III  in 
Figures  3a  and  3b, no  collisions  are  possible  since  the  surface  atom  is 
moving  away  from  the  gas  atom  at  t  ■  0,  On  the  plateau  and  for  most  of 
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the  apace  off  the  plateau,  region  I  In  Figures  3a  and  3b,  the  gas  atom 
strikes  the  surface,  returns  to  the  well,  and  rebounds  off  the  attrac¬ 
tive  step  before  another  collision  begins.  The  dashed  circle  in  Figure 
3b  indicates  the  veglon  of  interest  in  collision  initial  condition  space 
for  calculation  of  steady  state  adatom  energy  distributions  as  deter¬ 
mined  by  the  equilibrium  transition  rate  described  in  the  next  section. 

It  shows  that  only  collisions  in  region  I  affect  the  steady  state 
reduced  populations  and  rate  constants. 

In  region  II  on  Figures  3a  and  3b,  the  gas  atom  strikes  the  sur¬ 
face,  begins  and  ends  a  collision,  but  strikes  the  surface  again  before 
returning  to  the  well  and  rebounding  off  the  attractive  step.  To  assess 
whether  these  "double"  collisions  occur,  the  final  conditions  of  the 
normal  mode  solution  are  used  to  define  the  initial  conditions  of  new 
trajectories  for  the  time  t  >  t  .  The  gas  atom  is  in  free  flight,  and 
the  surface  atom  motion  is  approximated  by  simple  harmonic  oscillation 

a 

subject  to  initial  conditions  given  by  z(t  )  and  z(t  ).  Let  T  -  t  -  T  . 

c  c  c 

Then  the  gas  atom  and  surface  atom  equations  of  motion  are 

y(T)  -  y(T  )  +  y(x  )t 
c  c 

,  .  ,  %  •  (37) 

z(t)  ■  z(T  )  cos  T  +  z(x  )  sin  T  • 

v  C 

If  z(t)  ■  y(x)  for  any  T  >  0,  a  double  collision  occurs.  Region  II  in 
Figure  4  defines  the  values  of  y  and  V  for  which  double  collisions  occur. 
It  is  clear  that  double  collisions  are  not  important  for  most  values 
of  the  similarity  parameters  except  at  very  low  energies.  (This  energy 
dependence  derives  from  the  (y/e)  '  radius  of  the  region  of  interest  in 
^8»  space  defined  by  the  1/e  point  of  the  equilibrium  transition  rate, 
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R(e',  e) .  Double  collisions  may  be  important  for  cases  where  y  =  1, 

v 1,  and  5  is  small,  e.g.,  the  sublimation  and  condensation  of  rare 

(41) 

gases  on  rare  gases. 

Figure  5  shows  sample  gas  and  surface  atom  trajectories  illustrat¬ 
ing  the  various  types  of  collisions.  Nondime ns ional  displacements  are 
given  on  the  left  with  corresponding  velocities  on  the  right.  The  abscissa 
is  nondimens ional  with  respect  to  o>  Figure  5a  with  v  *  3.16  and  y  =  .5 

G 

represents  a  relatively  adiabatic  collision  at  Vc  *  1  and  Vg  =  0  in  region 
I  (single  collision)  on  Figure  3b.  Figure  5b  with  v*y*l,  V«0,  and 

V  ■  .2  is  on  the  boundary  between  region  I  (single  collision)  and  region  II 
(double  collision)  on  Figure  3a.  In  this  limiting  case  the  initial  condi¬ 
tions  for  Figure  5b  are  on  the  collision  duration  discontinuity  in  Vc, 

V  space.  Figure  5c,  with  v  »  y  «  1,  V  ■  0,  and  V  -  .86,  shows  a  rela- 

s  c  s 

tively  impulsive  double  collision  in  region  II  on  Figure  3a. 

IV.  GENERAL  SOLUTIONS 
A.  Equilibrium  Transition  Rate 

The  fundamental  role  of  the  one-way  equilibrium  transition  rate, 

R(e' ,  e)dede’  is  evident  from  the  functional  form  of  the  master  equation, 

Equation  (3),  and  the  adsorption  and  desorption  rate  constants  given  by 

Equation  (15)  and  (16)  where  R(f,  e)  is  defined  in  terms  of  R(e’,  e) 

by  Equation  5.  The  one-way  differential  transition  rate  may  be  expressed 

(38) 

as  a  flux  in  phase  space  from  the  state  e  to  the  state  e':v  ' 

dR  ■  p(v  •  n)ds  (38) 

where  ds  is  a  differential  element  of  a  surface  in  phase  space  separating  the 
states  e  and  s’  on  which  (v  •  n)  >0,  v  is  the  generalized  velocity  of  a  point 
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in  phase  space,  n  Is  the  unit  outward  normal  to  ds,  and  p(<^,  g)  is  the 

density  of  points  representing  the  state  of  a  system  of  n  particles 

in  a  6n  dimensional  phase  space,  the  axes  of  which  are  the  conjugate 

momentum,  ]>,  and  position,  £,  coordinates  of  the  particles. 

For  the  case  of  a  gas-atom,  surface-atom  collision,  the  only  non- 

lgnorable  coordinates  are  the  momentum  of  the  gas  atom,  p  ,  the  momentum 

g 

of  the  surface  atom  p  ,  and  the  position  of  the  surface  atom  Z.  The 

S 

one-way  equilibrium  transition  rate  may  therefore  be  written  in  terms 
of  the  initial  gas  atom  energy,  £,  and  the  initial  velocity  and  posi¬ 
tion  of  the  surface  atom,  V  and  V  ,  using  Equation  (38): 

C  8 

R(F’  Vr>  V  dFdv  dV«  "  P.»  ZX*  “  Z>dp  dp  dZ  .  (39) 

C  8  Co  js  B  fi  8 

The  surface  in  phase  space  across  which  the  flux  is  monitored  is  given 

•  • 

by  Y  -  Z  ■  0  with  the  condition  Y  -  Z  <  0. 

•  • 

With  the  definitions:  dp  ■  dE  /2u,  dp  ■  m  dZ,  and  Y  -  -  u, 

g  g  s  s  ’ 

Equation  (39)  becomes 

R(e,  Vc,  Vg)dedVcdVg  -  0(E?,  Z,  Z)  2ag(l  +  |)dEgdZdZ  .  (40) 

Since  it  is  assumed  that  the  gas  atom  and  surface  atom  are  Independent 
prior  to  collision,  the  phase  space  density  in  the  initial  state  p(E,  z,  z) 


can  be  separated  as  a  product  of  the  gas  particle  density,  o  (E  )  and 

P  g 

the  surface  particle  density,  p  (E  ): 

8  8 


P(E  ,  Z,  Z)  -  p  (E)  P  (E  )  - 
g  g  g  s  s 


N  exp  (-  E  /kT)  exp(-  E  /kT) 


,x - - -  (41) 

8 


where  N  is  the  number  of  gas  atoms,  Q  is  the  gas  partition  function, 
g  g 

Q  ■  2ricT/w  is  the  solid  atom  partition  function,  and  it  is  assumed 

8  S 


that  for  the  equilibrium  transition  rate  calculation,  the  phase  space 
points  have  a  Boltzman  distribution.  Substituting  Equation  (41)  into 
Equation  (AO)  with  Eo  =  kTe,  Ec  =  ^  mgu  (V^  +  Vr),  dZ  «  udV^,  and 


dZ  =  uu  ^"dV  yields 
s  s 


R(e,  V„,  VjdedV„dV„  ■>  ^  e"e  (-|— )  e“  u(Vc+V  (1  +  Vc)dedVcdVs  . 

4g 


c  s 


(A2) 


The  one-way  equilbrium  transition  rate,  nondimensionalized  with  respect 

to  N  kT/Q  is  then 
g  g 

^(e,  Vc,  Vg)dedVcdVs  =  d  +  Vc)  exp{-  ^  (V*  +  V*  +  tOJdedV^  (A3) 

where  the  nondimens ional  initial  gas  atom  kinetic  energy,  e,  surface 

atom  velocity,  V  ,  and  surface  atom  displacement,  V  ,  are  defined  in 
c  s 

the  previous  sections. 

from  the  definition  of  the  partition  function. 


NfkT  HK 

Qg  "  Qf  "  ~Q' 


(AA) 


where  the  subscript  g  indicates  all  gas  atoms,  the  subscript  f  indi¬ 
cates  only  free  atoms,  and  no  subscript  refers  to  adatoms.  Thus  the 
transition  rate  given  by  Equation  (A3)  may  be  nondimensionalized  in 
a  variety  of  ways.  In  adsorption  processes,  Equation  (A3)  is  nondi- 
mensional  with  respect  to 


NfkT 


j^.8kTv  1/2 

A'mn  J 
g 


-5 

e 


F  c 

&- 

A 


-5 

e 


—  •  1/2 

where  c  is  the  mean  speed  of  the  free  atoms,  and  Qf  ■  (27im^kT) 

for  gas  atoms  moving  normal  to  the  surface.  The  transition  rate 


(A5) 
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in  desorption  processes  is  nondimensionalized  with  respect  to 


_  A  co 

0  2tt 

where  Ag  is  the  equilibrium  adatom  concentration. 

Equation  (43)  shows  that  the  onc-’»ay  equilibrium  transition  rate 

has  a  strong  maximum  in  initial  condition  (V  ,  V  )  space.  In  the  limit 

c  s 

e  •*  V  ,  V  )  resembles  a  delta  function  at: 

V  =  0 


s 

Vc  =  j[(l  +  2p/e)1/2  -  1]  -►  0  .  (47) 

This  behavior  is  useful  in  making  the  variable  transformation  from 
/J(e»  Vc*  Vg)  to  J?(e',  e).  The  cnergv  transfer,  and  therefore  the  final 
energy  for  a  given  initial  energy,  is  given  as  an  implicit  but  exact 

function  of  e,  Vc,  and  Vg  by  Equation  (36).  This  function  is  plotted  in 

Figure  3b  for  the  particular  but  typical  case,  v  -  3.16  and  p  ■  .5.  It 
is  seen  that  the  surface,  ^  ( V c>  V  ),  is  approximately  a  plane  in  the 

neighborhood  of  the  maximum  of^(e,  V  ,  V  ).  The  equation  of  that  plane 

is  obtained  by  expanding  £  (V  V  )  in  a  two-dimensional  Taylor  series 
about  the  Vc,  Vg  values  given  by  Equation  (47)  and  truncating  the  expan¬ 
sion  at  the  linear  terms. 


e  <V  Vs*  "  a  "  1  +  b  vs  +  c  Vc  (48) 

where  a  -  y2(tc),  b  =  2y(Tc)  yy  (tc),  and  c  -  2y(tc)  yy  (t  )  with  sub- 

s  c  c 

scripts  indicating  differentiation,  and  Tc(0,  0)  is  obtained  from  Equa¬ 
tion  (33).  A  rotation  transformation  to  two  new  coordinates,  V  and 

’  sr 

vcr»  whose  axes  are  parallel  and  normal  to  the  intersection  of  the  e' 
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plane  with  the  plane  of  the  V  ,  V  axes,  picks  up  the  direction  of  the 

C  8 


normal  coordinate  V  .  This  transformation  is 

cr 


-  vad  +  V 

c 


(49) 


2  -1/2  .  v2  -1/2 

V  -  V  (1  +  £*) 
sr  c  jT 

.2  -1/2  2  -1/2 
Vcr-Vl+V  +V.<1  +  J>  • 

c  b 

The  Jacobian  is  unity.  The  V  ,  V  axes  can  be  located  on  Figure  3b 

cr  sr 

by  an  ^  60-degree  counter-clockwise  rotation  from  the  V£,  axes. 
Inverting  Equations  (49),  substituting  into  Equation  (43),  and  integrat¬ 
ing  over  V  gives 

S  t 

1/2  2  -1/2 

R  <V  E>  ■  <5T>  «  +  'cr'1  +  V  1  “»<-  ?  Vcr  -  E'dVcrdE  ' 

(50) 


The  one-way  equilibrium  transition  rate  is  then  obtained  from  Equation 

I  dV 

(50)  using  K  (e,  e')  de'de  -  /i(V  (e'),  e) 


cr 


cr 


de' 


de'de  and 


Vcr(e',  e)  ■  (e'  -ae)e~*  (b2  +  c2)”1^2  with  the  desired  result 

P  tr  ri^  .  e(c2  ±  b2)  +  c(e'  -  ae)  ..._f-(e'  -  ae)  >  rsn 

A  U,  C  )  ■  .  2  2  37?  2  i  “  (51) 

(ny)1'2  (e (c  +  bz)r'z  ye(bz  +  cz) 

This  approximate  rate  does  not  satisfy  the  symmetry  required  by  detailed 

balancing.  K  symmetric  form  can  be  constructed  by  taking  the  arithmetic 

mean  of  R(c,  e')  and  R(c',  e).  Numerical  comparisons  show  negligible 

distinction  between  the  symmetric  and  unsymmetrlc  rates. 

Introducing  the  approximations:  a  ^  1,  e  e  ■  (e  +  e')/2,  and 

A/2e  «  1  yields  a  simple  symmetric  kernel: 


R  (e,  A)  - — T7T  exp  ^ - 1~2 - T~  ~  ^  (52) 

(ywe(bz  +  c))1'2  ye(b  +  c  ) 

where  e  •  (e'  +  e)/2  is  the  mean  energy,  and  A  -  c'  -  c  is  the  energy  trans¬ 
fer.  This  approximate  one-way  equilibrium  transition  rate  is  used  in  the 
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calculatlons  reported  here.  Figure  6  shows  contours  of  constant  non- 
dimensional  one-way  equilibrium  transition  rate c')  in  the 
range  10  to  10  in  c,c'  space.  The  gradient  of  /i(c,  e')  along  the 
A  axis  is  much  greater  than  along  the  c  axis.  This  difference  allows 
an  expansion  of  the  population  X(c,  t)  about  A  ■  0  to  obtain  an 
equivalent  diffusion  equation  for  X(c,  t)  frem  the  master  equation. 

The  exponential  decay  in  c  follows  from  the  Boltznan  factor.  A  typi¬ 
cal  well  depth,  6  •  13,  is  indicated  by  the  dashed  line  in  Figure  6. 
There  are  four  transition  regions  separated  by  c*  •  6  and  c  •  6:  for 
c  and  e'  >  6,  the  gas  atom  is  never  trapped;  when  C  and  c*  <  6,  transi¬ 
tions  occur  between  states  in  the  well;  for  C  <  6  and  s*  >  6  desorption 
takes  plsce;  and  in  the  region  C'  <  6  and  e  >  6  gas  phase  atoms  are 
adsorbed. 

Once  the  equilibrium  transition  rate  is  km*wn,  the  one-way  equi¬ 
librium  desorption  rate,&(f,  c),  is  obtained  as  the  Integral  of 
JKe',  c)  over  6  <  e'  <  "  given  by  Equation  (5).  This  is  equivalent  to 
integrating  Equation  (52)  over  .5(6  +  c)  <c<«*  using  A  -  2(c  -  c) 
wit h  c  as  a  parameter.  The  result  is 

*(f,  e)  «  -  exp(-  4  qpc  +  8q2c)  lerfc(pi  -  2|^)  +  e8<,pC  erfc(pi  +  -|£)1 

(53) 

where  q  -  (p(b2  ♦  c2))’1/2,  p  -  (4q2  +  1)1/2,  t  -(~^)1/2  . 

The  Causslan  behavior  of  the  kernel  for  Equation  (5)  given  by 
Equation  (52)  is  clear  in  Figure  7  where  the  nondlmensional  transi¬ 
tion  rate  is  plotted  versus  energy  transfer  for  w  •  .5  and  v  -  3.16. 

The  two  curves,  equivalent  to  slicing  Figure  6  along  c  ■  10  and  c  -  30, 
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lndicatc  thsc  the  cean  energy  transfer  Increases  as  che  Incident  energy 
Increases.  This  Is  reasonable  since  harder  collisions  should  be  sore 
impulsive.  The  broadening  of  che  kernel  half-vldth  appears  In  Equa¬ 
tion  (52)  as  a  Causslan  scale  factor  of  (c)”*. 

8.  Steady  State  Adacoa  Energy  Distribution 
Two  methods  were  used  to  obtain  approximate  solutions  f,e  the  caster 
equation  for  the  steady  state  adatoa  energy  distribution. 

t.  Equivalent  Diffusion  Equation 

Equation  3  may  be  converted  to  an  equivalent  diffusion  equation^7) 

which  identifies  the  desorption  process  with  a  diffusion  process  In 
(24) 

phase  space.  This  is  a  particular  form  of  the  sore  general  crans- 

(25) 

formation  to  a  Fokkr  r-Pianck  equation.  In  che  steady  state  approxima¬ 

tion  the  equivalent  diffusion  equation  is 

§f  (~2 - «  o  (54) 

where  A  (c)  *  (e*  -  e)n/J (c,  c')  de*  is  the  moment  of  tie  energy 

transfer  with  respect  to  the  equilibrium  transition  rate.  Equation  (54) 
is  a  valid  approximation  to  the  master  equation  when  the  form  of/?(c,  c') 
is  such  that  ®  j  3A2/3c. 

Figure  8,  where  v  •  3,16  and  y  •  .5,  indicates  this  criteria  is 
satisfied  for  c  >  2  fot  the  kernel  given  by  Equation  (52).  This  result 
is  typical  for  all  values  of  the  Inertia  and  frequency  parameters.  Fig¬ 
ure  8  also  displays  the  exponential  energy  dependence  of  Aq  (c)  due  to 
the  Boltzmnm  factor  ir/£(C,  c').  The  normalization  condition,  x(0)  *  1# 
replaces  the  boundary  condition  at  Che  bottom  of  the  well  since  the  solu¬ 
tion  of  Equation  (54)  is  not  well  defined  with  first  derivative  boundary 
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A«  W)  j/iv 

conditio  a*.  Th«  boundary  condition  it  c  ■  6  it  — j •  -  2T  (6)  xW) 
where  6  le  the  well  depth  end  T  (6 )  it  the  "one-way  equilibrium  desorp¬ 
tion  flux": 

r(6)  -  /£/?(c,  c*)  dc’  dc  .  (55) 

Once  the  equilibrium  transition  rate  is  known,  T(6)  can  be  integrated 
directly.  Using  the  approximate  kernel  given  by  Equation  (52)  and 
changing  variables  to  c  and  6,  T(6)  is  obtained  by  integrating  over 
2(6  -  c)  <  A  <  «*  for  -  •  <  c  <  6  aqd  over  2(c  -  6)  <  A  <  •  for  6  <  c  <  «* 
with  the  result 

T(6)  -  exp (-6  ♦  6/16q2)  erf(61/2/4q)  .  (56) 

26 

This  approximate  equilibrium  desorption  flux  is  valid  when  a  «  1,  ~  »  1 
and  8q2  »  1. 

The  steady  state  reduced  population  in  the  diffusion  approximation 
is,  integrating  Equation  (54)  subject  to  the  above  boundary  conditions: 

X(C)  •  1  -  /J  A2(g)"1  dc/l/jA^c)**1  dc  +  \  r(6)“lJ  .  (57) 

Numerical  integration,  uelng  Equation  (56)  for  T(6)  and  A^fc)  as  given 
in  the  Appendix  gives  the  results  described  here. 

The  diffusion  approximation  is  useful  since  its  range  of  validity 
Includes  significant  nonequlllbrlua  effects,  and  the  moments  of  the 
transition  r av«*#?(c,  V  ,  V  )  as  shown  in  the  Appendix  are  sufficient 
to  obtain  a  solution.  A  criterion  for  validity  of  the  diffusion  solu¬ 
tion  is  given  by  the  mean  square  energy  transfer  per  collision  A2(c)/AQ(e), 
shown  in  figure  9  for  v  •  3.18,  y  «  .72,  and  v  •  3.14,  y  »  .46.  When 
A9/Aq  <  1,  the  diffusion  solution  is  accurate.  When  A2/Aq  >  1*  che 
teres  0(A)  dropped  In  the  diffj’  jion  approximation  become  Important. 
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Flgurc  10a  shows  equivalent  diffusion  equation  solutions  given  by 
Equation  (57)  for  the  steady  state  reduced  population  for  v  ■  3.16, 

U  •  .5,  and  6  «  10  and  30.  The  ratio  of  the  actual  to  the  equilibrium 
energy  distribution  is  plotted  versus  energy  below  the  top  of  the 
well.  Those  energy  levels  for  which  x(c)  equals  one  have  a  Boltmann 
distribution.  The  levels  for  which  X(£)  is  less  than  one  are  depleted. 
The  magnitude  of  the  depletion  is  determined  by  the  rate  at  which  atoms 
"diffuse"  up  the  well  to  replace  the  adatoms  that  desorb.  When  the 
energy  transfer  in  a  collision  is  small,  diffusive  movement  of  adatoms 
among  the  energy  states  in  the  well  is  slower  than  the  desorption 
process, and  the  upper  levels  deplete.  However ,  when  the  energy  trans¬ 
fer  is  large,  the  desorbed  atoms  are  rapidly  replaced  by  adatoms  from 
lower  levels,  and  the  true  population  approaches  an  equilibrium  popula¬ 
tion. 

2.  Iteration  Solution 

When  the  energy  transfer  per  collision  is  large,  A^/Aq  >  1.  itera¬ 
tion  schemes  converge  rapidly,  and  the  first  iteration  is  a  valid 
approximation  to  X(c).  The  steady  state  adatom  energy  distribution 
is  then  obtained  from  an  iteration  solution  to  Equation  (10)  in  the 
form 

Xj(e)  -  /J*(e,  g,)Xj_1(c')  deV/JJ*  (c,c')  dc'  (58) 

t»h 

where  x^ (e)  is  the  J  iteration  approximation  to  the  steady  state 

reduced  population.  Assuming  the  trial  function  Xq(c)  “1.0  for 

C  <  6  and  substituting  Equation  (52)  into  Equation  (58)  with  a  change 

-  1/2 

in  integration  variable  from  c'  to  (c)  gives 


(59) 


Xi<e)  •  /J*  (*,  c)  dx/y£(R(x,  c)  dx 

2  7  2  2  —2  «•  1/7 

where  /((x,  C)  -  —  exp  -{p  x  -  R  x  +  pr)  with  x  •  (e)  '  and  r  •  2qc. 

/? 

Integration^**^'  yields  a  first  lnteratlon  approximation  for  the 
steady  state  reduced  population. 

Xjfc)  -  1  -  >(erfc(pt  -  ^)  +  erfc(pt  +  exp(Bpqc))  .(60) 

A  second  Iteration  approximation  for  the  steady  state  population  at  the 
top  of  the  well,  x2(^*  »ay  be  found  by  msserical  Integration  of  Equa- 
tion  (58)  with  J  •  2  and  C  •  6  using  Equation  (59)  for  Xj.jtc').  The 
analytic  single  iteration  approximation,  Equation  (60),  is  used  in 
calculations  reported  here.  Although  no  direct  comparison  is  possible 
Figure  10b  shows  the  equivalent  diffusion  equation  and  iteration  solu¬ 
tion  approximate  steady  state  reduced  populations  for  v  *  3.18,  U  •  .72 
at  6  •  55.6.  X  indicates  a  second  iteration,  x2(6).  <n*e  "*a!X  energy 
transfer  per  collision,  Aj(6)/Aq(6),  is  1,5.  The  actual  population 
is  close  to  equilibrium,  and  similar  ratios  of  the  steady  state  to 
equilibrium  rate  constants  are  obtained  from  both  methods;  Diffusion 
equation  *  .86,  and  Iteration  solution  kj/hjg  *  *88. 

€,  Steady  State  Pate  Constants 

The  energy  Integral  of  the  master  equation  yields  the  adsorption 
and  desorption  rate  constants  in  the  steady  state  approximation.  Two 
methods  have  been  described  for  obtaining  the  steady  state  reduced 
population  and  the  equilibrium  desorption  flux,  required  bv  Equations 
(15)  and  (16)  for  these  rates.  The  iteration  solution  msserically 
integrates  Eouations  (15)  and  (16)  using  Equation  (60)  for  the  stcadv 
reduced  population,  X(E)*  and  Equation  (53)  for  the  equilibrium  desorption 
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race,  ^(f,  c).  The  diffusion  method  formally  integrates  Equations  (IS) 
and  (16)  t£ing  Equation  (57)  for  the  steady  state  reduced  population, 
X(c),  and  Equation  (56)  for  the  desorption  flux,  T(6),  required  as 
a  bounda-y  condition.  Since  the  iteration  method  la  valid  over  a 
restricted  range  of  similarity  parameters  where  the  adatom  enacefclc 
is  close  to  equilibrium,  the  diffusion  meth\>d  is  used  for  the  results 
described  here. 

The  steady  state  desorption  rate  constant,  obtained  in  the  diffu¬ 
sion  approximation  froo  Equation  (15)  noting  that  NkT/Q  »  A  u  /2n  for 

e  g 

the  collision  model  described  in  Section  III  is 

kd  *  l/o  a2X  dc  ♦  ?  rw)*1))'1  (61) 

A 

where  Aj  is  given  by  Equation  (A9)  and  T(6)  is  given  by  Equation  (57). 
The  equilibrium  desorption  rate  constant  is  obtained  from  Equation  (15) 
using  x(c)  •  ui 


(62) 


The  classical  desorption  rate  constant,  kj^,  is  given  by  Equation  (2). 

The  classical  derivation1  assumes  that  (1)  the  gas  phase  and  the  adsorbed 
phase  are  at  equilibrium;  and  (2)  the  sticking  coefficient  is  unity  so 
chat  the  rate  of  arrival  equals  the  rate  of  adsorption.  These  assump¬ 
tions  are  not  physically  realistic  during  the  adsorption-desorption 
relaxation  processes.  The  equilibrium  desorption  rate  constant.  Equa¬ 
tion  (62),  is  valid  when  the  first  assumption  is  correct,  and  the 
"steady  state"  desorption  rate  constant.  Equation  (61),  does  not  require 


either  assumption.  For  reference  the  equilibrium  constant  for  the 
adsorption-desorption  process  is 
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K  •  (■  /2nkT)1^2  «  exp  (-  6),  (63) 

*  &  f 

end  the  classical  adsorption  rate  constant  is  the  arrival  rate, 

k  -  -  (kT/2wal1/2  .  (64) 

ac  4  *  g 

When  Equation  (16)  la  valid,  the  equllibrlua  adsorption  rate  constant 
la 

k  -  (kT/2im  J1/2  T(6)  exp  (6)  ,  (65) 

s 

and  the  corresponding  steady  state  adsorption  rate  constant  la 

ka  -  lkT/2ir»g)1/2  (2  ^(c)"1  dc  +  \  rW)-1)'1  exp  (6). 

g  v  (66) 

Equation  (17)  shows  that  the  ratio  of  steady  state  to  equilibrium 
rate  constants  Is  Independent  of  the  direction  of  the  process.  Froa 
Equation  (17)  with  Equations  (61)  and  (62)  or  Equations  (66)  and  (65), 
k  k 

\  +  21* (6)  /*  AAcf1  dc)"1  •  (67) 

ae  de 

k  kd 

In  the  adiabatic  Halt  A2  *♦  0,  and  - —  ■  r — ■*  0.  In  the  lepulalve 

ae  ae 

Halt  A2  la  a  aaxlsua.  Since  v  «  gives  the  adiabatic  Halt  and 
6  <*  gives  the  lapulalve  Halt,  nonequlllbrlua  effects  generally 

increase  as  V  Increases  and  decreases  as  6  Increases.  Figure  11 
presents  the  ratio  of  the  steady  state  tc  equlllbrlua  rate  constants 
given  by  Equation  (67)  for  u  «  .2  and  1  with  6  •  10.  Under  this  condi¬ 
tion  of  fixed  6,  the  classical  desorption  and  adsorption  rate  constants 
for  a  given  gas  specie  would  be  constants. 

An  interesting  effect  Is  shown  by  the  contours  of  constant  u  in 
Figure  11  which  pinpoint  the  values  of  the  siallarity  paraseters 
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corrcspondlng  to  highly  nonequilibrium  adatom  energy  distributions. 

"Resonance"  occurs  when  the  larger  normal  mode  frequency  given  by 

Equation  (30)  is  an  odd  multiple  of  the  smaller: 

Uj  •  (2n  +  Du^  :  n  -  1,  2,  3  ....  (68) 

For  normal  mode  frequencies  which  are  not  independent,  the  period 

of  both  the  gas  and  the  surface  atom  motion  is  a  multiple  of  the 

longer  of  the  two  normal  mode  ocrlods,  The  end  of  a  collision 

occurs  at  t  ■  t  ■  ir/u^,  and  since  y  and  z  have  a  co^mn  period,  the 
•  • 

velocities  y(Tc)  and  z(tc)  return  to  their  initial  magnitudes.  Thus 
at  resonance  the  energy  transfer  is 

A/c  -  (y  (t£) 2  -  y(0)2)  -  0  ,  (69) 

and  the  adiabatic  limit  is  reached.  At  frequency  parameters  near 
resonance, the  steady  state  rate  constant  is  orders  of  magnitude  smaller 
than  the  equilibrium  rate  constant.  Resonance  is  therefore  a  special 
case  where  non^qul librium  effects  are  extremely  important.  Including 
enharmonic  damping  would  round  off  the  siinlma  in  Figure  11,  but  would 
not  change  the  qualitative  picture. 

Using  Equation  (30)  for  the  normal  mode  frequencies  in  Equation  (68) 
obtains  the  resonance  values  of  v  as  a  function  of  u  plotted  in  Figure 
12; 


vr0i)  -  (2A’(1  ♦  (1  -  $(  1  ♦  p))1'2)  -  (1  +  u)) 


where 


*  .  1  .  2.  , 

( (2n  +  1  r  +  1) 


,1/2 


(70) 


Vf(u)  goes  to  2n  +  1  in  both  the  limit  n  ",  and  the  limit  li  0.  With 
these  results  and  the  steady  state  adatom  energy  distribution  previously 
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derlved,  the  magnitude  and  regions  of  import  of  nonequllibrlua  effects 
have  bsen  defined. 

^Iure  shows  the  variation  of  the  adsorption  rate  constants 
given  by  Equations  (65)  and  (66)  for  v  -  3.18,  y  •  .72,  tad  v  *  2.66, 
y  -  .38  over  a  range  of  distension  less  well  depths  6  fron  10  to  80. 

The  adsorption  rate  constant  provides  a  aore  accurate  measure  of  non¬ 
equllibriua  sffects  sines  the  exponential  teaperaturc  dependence  of 
the  desorption  rate  constant  is  not  present.  The  steady  state  and 
equilibrium  rate  constants  shown  are  nondiaenslonallred  with  respect 
to  the  arrival  rate,  and  the  nondiaenslonal  well  depth,  6,  is  equivalent 
to  an  Inverse  teaperaturc.  The  dlsplaceaent  of  the  k  lines  below 

AG 

kac  *  1  *iv*»  a  aeasure  of  the  effect  of  a  non-unity  sticking  coeffi¬ 
cient.  The  corresponding  dlsplaceaent  of  the  k  lines  below  the  k 

®  AG 

lines  indicates  the  effect  of  a  nonsqulllbrlta  surface  population. 

The  results  given  here  are  for  adsorption  fron  an  Incident  gas 
with  a  Bolt i*an  energy  distribution  about  the  surfacs  teaperaturc. 

When  the  incident  gas  streaa  has  any  other  energy  distribution,  the 
adsorption  rate  constant  is  obtained  fron  Equation  (19).  The  result 
is  Equation  (66)  multiplied  by  H(0)  where  H(0)  is  given  by  Equation 
(9)  as 

H(0)  •  <^/?(C,  0)  X£C)  de//£/?  (c,  0)dc  .  (71) 

for  the  case  of  a  gas  with  a  Boltnean  energy  distribution  at  a  tcaperi.- 
1/2 

ture  V  *  n  exp*~  H(C  -  6)  +  c  -  6)  where  n  -  T/T  .  Substituting 
this  expression  and  the  value  of  K  (c,  0)  obtained  from  Equation  (51) 
into  Equation  (71)  and  Integrating  gives 
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h(o^  -  n  '  texp«(n  -  D)  - > 


V  +  Ul/2  erfc[6(a2g2  +  n))1/2 


a  q  +  n 


erfc[6(a  q  +  1)] 


For  Che  case  of  a  oonocnergetlc  incident  beaa  of  energy  Cq, 

)^(c)  -  (4VQ/c)  6 (c*6~Gq)  exp(e  -  6),  and  Equation  (71)  gives 

.....  %  r.V  +  i  ,i/:  cxp['  ‘^o  -  <*2"2  +  1)41  .... 

H<°>  *  ~  <w?i-r5T>  — TSSuTTIw1"-  (73> 

1/2 

whcr<*  V„  •  (2kTc_/m  )  la  the  velocity  of  the  atoms  in  the  bean, 
u  u  g 

_  1/2 

and  c  -  (SkT/iio  )  in  the  Man  themal  speed  for  a  gas  of  a  tots*  in 
g 

equilibrium  with  the  surface.  Figure  14  shows  H(0)  n  *^2  given  by 
Equation  (72)  versus  the  frequency  parameter,  v,  for  n  ■  0,  1,  and 
10  with  y  •  .2  and  1.  The  wall  depth,  6,  was  fixed  at  d  ■  10;  there 
was  no  Masurable  variation  of  H(0)  with  well  depth  over  10  <  6  <  30. 
H(0)  o"1^2  •  1  for  all  v^(y)  and  for  n  •  1. 

-1/2 

figure  15  shows  the  n  dependence  of  11(0)  n  ftiven  by  Equa¬ 
tion  (72)  for  '■  •  3.18,  y  «  .72,  and  v  ■  2.66,  y  •  .38.  The  effect 

of  tho  former  set  of  paraMters  being  "off  resonance"  while  the  latter 

-1/2 

is  "on  resonance"  is  clear  from  the  near  unity  H(0)  n  for  the 


resonance  system. 

As  a  sample  calculation  of  the  steady  state  rate  constants, 
consider  the  case  of  a  xenon  gas  impinging  on  a  clean  tungsten  sur¬ 
face.  Comparison  of  the  soft  cube  model  with  scattering  pattern 
(18) 

data  Indicates  a  frequency  parameter  v  •  2,0  and  a  well  depth 
D  •  4500  °K  for  this  system.  The  inertia  ratio  is  U  m  .72.  If  the 
surface  temperature  were  100  °K,  6  would  be  45,  If  T  ■  100  °K, 

3  -1 

n  •  1,  and  H(0)  »  1;  therefore,  from  Equation  (66)  kfl  »  2.0  x  10  cm  sec  . 
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If  Tg  -  1000  °K,  n  -  0.1  and  from  Equation  (72)  H(0)  -  .36;  therefore, 

from  Eq  in t ion  (19)  kQ  -  .72  x  103  cm  sec"1.  The  desorption  rate 

constant  is  the  isamc  for  all  T  ;  here  k.  -  4.5  x  10“8  sec’1.  Note 

8  ® 

that  the  energy  distribution  of  the  incident  stream  primarily  affects 
the  sticking  coefficient  while  the  steady  state  distribution  given 
by  Equation  (57)  remains  a  good  approximation  to  the  adatom  energy 
distribution  for  all  r\. 

V.  CONCLUDING  REMARKS 

In  summary,  the  time  evolution  of  a  nonequilibriua  ensemble  of 
gas  atoms  adsorbed  on  a  solid  surface  has  been  described.  The  tem¬ 
perature  and  atomic  species  dependence  of  the  phenomenological  rate 
constants  in  the  experimentally  observed  first-order  rate  equation 
has  been  obtained  from  the  energy  integral  of  the  governing  master 
equation.  The  pertinent  similarity  parameters  and  energy  transition 
probabilities  for  gas  atoms  colliding  with  solid  surfaces  have  been 
derived  from  a  classical  interaction  (model  using  a  gas-surface  poten¬ 
tial  composed  of  an  oscillating  harmonic  repulsion  and  a  stationary 
attraction.  The  master  equation  for  the  reduced  adatom  energy  popula¬ 
tion  has  been  solved  in  the  steady  state  approximation  by  intcration 
when  the  mean  energy  transfer  per  collision  is  greater  than  kT  and 
by  conversion  to  an  equivalent  diffusion  equation  with  the  second 
moment  of  the  energy  transfer  assigned  the  role  of  diffusion  coeffi¬ 
cient  in  the  case  of  mean  energy  transfer  less  than  kT. 

The  magnitude  of  the  departure  from  equilibrium  due  to  the  deple¬ 
tion  of  the  energy  levels  near  the  top  of  the  surface  potential  well 
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is  defined  by  the  steady  state  reduced  adatom  energy  population.  It 
wss  found  that  for  "resonance"  values  of  the  ratio  of  surface-to-gas 
oscillator  frequencies,  gas-surface  collisions  were  sufficiently 
adiabatic  that  the  adatoa  energy  distributions  were  sppreclably  non- 
equlllbrlua,  end  the  stesdy  state  rate  constants  were  depressed  orders 
of  magnitude  below  the  equlllbrlua  rste  constants.  Experloental  verifica¬ 
tion  of  this  prediction  would  be  extrenely  useful. 

The  ratio  of  the  stesdy  state  desorption  to  adsorption  rate 
constants  is  the  equilibrium  constant,  unless  an  incident  gas  exists 
with  an  energy  distribution  which  is  not  Boltzman  about  the  surface 
teaparaturc.  If  such  s  gas  exists  the  adsorption  rate  constant  is 
modified  by  the  ratio  of  the  adsorption  rate  at  the  bottom  of  the  gas- 
surface  potential  well  from  the  actual  incident  energy  distribution 
to  the  adsorption  rate  at  that  energy  level  from  an  equilibrium  inci¬ 
dent  energy  dlatrlbutlon.  The  desorption  rate  constant  is  unmodified; 
therefore,  the  ratio  of  rate  constants  is  the  product  of  the  equilibrium 
constant  and  the  ratio  of  adsorption  rates  at  the  bottom  of  the  well. 

Comparisons  with  experimental  data,  as  described  in  a  later 
paper,  indicate  the  results  of  this  analysis  are  reasonable;  however, 
insufficient  experimental  results  are  available  at  the  present  time 
to  say  that  the  theory  has  been  confirmed.  This  treatment  obtains 
the  detailed  time  history  of  the  energy  distribution  of  the  adsorbed 
ensemble.  The  adsorption-desorption  rate  problem  considered  here  is 
but  one  of  a  large  set  of  problems  whose  solution  may  be  described  by 
similar  analysis.  Gas-surface  potential  parameters  obtained  from  this 
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approach  penal t  prediction  of  many  phenomena;  therefore,  rigorous 
experimental  verification  la  feasible. 

Possible  extensions  of  the  model  are  (1)  to  account  for  surface 
roughness  by  introducing  a  distribution  of  surface  orientations  or 
well  depths;  (2)  to  Include  energy  transitions  Induced  by  adatom-adatom 
collisions  in  the  master  equation,  thus  eliminating  the  low  surface 
concentration  restriction;  and  (3)  to  add  another  degree  of  freedom 
to  the  master  equation  to  account  for  internal  energy  exltatlon,  thus 
extending  the  analysis  to  adaolecules. 
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APPENDIX:  ENERGY  TRANSFER  HOMENTS 


The  Doacnts  defined  In  Section  IV  ere  evaluated  here  for  the 

collision  model  described  In  Section  III.  Using  the  equilibrium 

transition  rate  given  by  Equation  (43*  for  the  required  kernel  and 

recognising  that  Integration  over  c*  Is  equivalent  to  integration 

over  V  and  V  ylel * 
c  % 

A  (c)  -  r  r  An  •  (1  +  V  )  exp  -  £  (V2  +  V2  +  u)dV  dV 

n  -*»  Tip  c  r  p  c  s  M  c  s 

(Al) 

where  the  lower  Unit  of  Vc  has  beji  extended  to  — »  with  a  negligible 
additional  contribution  to  the  Integral.  Since  li;:c  energy  transfer , 

A.  Is  implicitly  dependent  on  \'c  and  V#  through  the  collision  deration 
Tc»  It  is  necessary  to  approxinate  A  by  a  two-dinenslonal  Taylor  expan¬ 
sion  about  the  naxlnua  In  the  kernel  taken  here  to  second  order.  This 
expansion  Is 

A(VcVs)  •  A  +  Ay  <Vc  “  Vc)  ♦  Ay  (V#  -  V#) 

C  8 

4'*vt  <*c-v2  +  “»v  <vc-v(v. -v 

c  c  vs 


+  Ay  y  (VB  -  v/j  (A2) 

where  the  tilde  superscript  Indicates  evaluation  at  the  kernel  maximum 

~  I/O 

(V  •  0,  V  ■  ,5((1  +  2p/e)  '  -1))  and  subscripts  indicate  differentia- 

8  C 

tlon.  Third-order  terms  arc  odd  and  contribute  negligible  to  the 
Integral.  The  coefficients  In  Equation  (A2)  are 


A 


* 


e(y2  • 

*  U  . 

Ay  y  -  2  e(yj 

+  y  y v  v  ) 

c  c  c 

C  C 

2  c  y 

• 

Ty  , 

• 

Ay  y  -  2  G(yv 

•  •  • 

7y  +  7  7V  v  ) 

c 

C  S  8 

c  c  s 

• 

2  e  y 

7y  . 

~  *2 
Ay  y  -  2  e(yy 

+  y  y v  ) 

8  8  8  8  8 


(A3) 
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where 

with 

y  •  Vl  +  Vcf2  +  f3  •  *V  ■  f2  •  •  fl 

c  s 

(A4) 

e 

^  W1  *ln  Vc  "  “2  *in  Vc^1  " 

*  2  2  2 
f2  •  (cos  U2Tc  -  COS  UjTc>  /  (u)j  -  Ujjv 

(A5) 

f3  •  ((1  -  V2Uj)  COS  -  (1  -  V2U2)  COS  U1T<;)/(w2  - 

w2)v2 

and  where 

yv  v  •  1(U1  *ln  Vc  “  w2  8in  Ve)/(ul  •  w2)N,2jTcV 

c  c  c 

yv  V  "  l(wl  *ln  Vc  “  U2  ,in  Vc)/<W1  “  u2)v2,TcV 
yV  V  “  lw2 (cos  MXT  -  COS  '*>2 0/(1  -  C) It  v 

S  S  B 

(A6) 

with 

TcV  "  <ul  *ln  Vc  •  u2  ,ln  Vc)/<wl  "  u2)(y  “  t)v2) 

(A7) 

TcV  "  w2(co*  uiTc  *  C0B  w2Tc>/((1  *  t)(y  "  r)) 
vhtr#  Tc(V#t  V£)  Is  evaluated  via  Equation  (33).  ‘these  expressions 


were  obtained  by  differentiating  Equation  (31). 

Substituting  Equa- 

tlons  (A2)  through  (A7)  Into  (Al)  and  Integrating  obtains  the  follow- 

ing  SOMDtS. 

Aq(c)  •  exp  (r  C) 

(A8) 

AX(C)  •  (exp  (•  €)){A  ♦  (i~  -  vfi)  ^ 

c 

(A9) 

c  c 

♦  ^  y  )) 

s  s 
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A2(c)  •  (e*p  (•  c>)  (A2  +  (jj  -  V£)  AAy  (A10) 

♦  h  ,(T  Vl  '  »c  *  11  y  I)  • 

c  C  C  S  SB 

The  zeroth  oocent  is  singly  the  collision  rate.  The  first  eoaent,  as 
indicated  in  Figure  8  satisfies  the  criteria*27*  Aj(c)  *  •  *  3A2*c*^8c 
required  for  validity  of  the  diffusion  approximation.  A  measure  of 
the  mean  energy  transfer  per  collision  is  given  by  the  ratio  of  the 
second  to  zeroth  moments.  The  doainance  of  second-order  teres  in  the 
energy  transfer  Taylor  expansion  is  exhibited  in  the  linear  dependence 
°f  A2^A0  on  C  a®  ®hown  ln  9.  Using  Equation  (A10)  to  evaluate 

in  Equation  (30)  gives  the  desired  steady  state  solution  to  the 
approximate  diffusion  equation  for  the  adatou  energy  distribution. 

In  addition  to  these  moments  an  approximate  one-way  equilibrium 
desorption  flux  is  given  by  the  first  moment  of  the  absolute  value  of 
the  energy  transfer 

T(c)  \  |fi|L  (e)  -  \rM  |A|  R(e,  e')de'  .  (All) 

Using  a  rotation  transformation  similar  to  that  used  to  obtain i\(.e,  e’) 
from  J?(e,V  V  ) ,  Equation  (All)  may  bo  separated  into  a  summation  of 
integrals  in  which  A  has  a  constant  sign,  thus  bringing  the  absolute 
value  outside  the  individual  integrals.  Asymptotically  expanding  the 

1 

resulting  parabolic  cylinder  functions  obtains 

\L\L(c)  -  .5e  (Aj  +  Aj  )1/2  exp (-  e  -  l|g£)  sinh  [(5e/y)1/2  C] 
sc 

2 

+  2  (1  -  BC)  (m/tte)1/2  exp  (-  cosh  [  (3e/u)1/2C]  } 


(A12) 
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where 


B  -  (1  +  aJ  /aJ  )*1/2 

8  C 

C  -  (A  -  V.Ay  )/(Aj  +  Aj 


1/2 
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Fig-  1 

Fig-  2 

Fig.  3 

Fig.  4 

Fig.  5 


FI CURE  CAPTIONS 

Gas-surfa.re  potential  veil  formed  by  superposition  of 
dynamic  repulsive  and  static  attractive  potentials. 

Collision  model.  Attractive  potential  Is  an  arbitrary  shape 
of  height  6  at  fixed  position  several  lattic  spacings  from 
the  surface.  Repulsive  quadratic  potential  follows  the 
surface  oscillator  motion. 

Contours  in  initial  condition  space:  (a)  contours  of 

constant  collision  times  for  v  ■  y  ■*  1;  (b)  contours  of 

constant  energy  transfer  for  v  -  3.16  and  y  «  .5.  Region  I 

has  single  collisions,  and  Region  II  has  second  collisions. 

In  Region  III  no  collisions  occur.  Heavy  dash  indicates 

plateau  bounded  by  a  discontinuity  in  t  and  A. 

Second  collisions  are  not  significant  at  the  incident  gas 

atom  energies  Indicated  when  the  system  frequency  and 

inertia  parameters  lie  in  Region  I  on  this  figure. 

Gasy(t)  -  ,  and  surface  z(t) - atom  trajectories. 

(a)  adiabatic  collision  (v  ■  3.16,  y  =  .5,  V  ■  0,  V  ■  1.0) 

s  c 

with  Tc  =  3.06  and  A/e  =  -  .007;  (b)  limiting  case  between 

single  and  second  collisions  (v  *  y  =  1,  V  =  .2,  V  =0) 

s  c 

with  t  =  2.80  and  A/e  *  -  1  at  the  end  of  the  first  colli- 
c 

sion;  (c)  second  collision  (v  =  y  =  1,  V  =  .86,  V  =0) 

s  c 

here  t  =1.71  and  A/e  =  -  .69  at  the  end  of  the  first  colli- 
c 

sion.  Length  and  time  are  nondimensional  with  respect  to 
u  wg1  and  ujg1,  respectively. 
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.  6  Contours  of  constant  one-way  equilibriums  transition  rate  in 
c,  c'  space  for  v  -  3.16  and  p  -  .5.  Saeple  well  depth 
shown  at  6  ■  13. 

.  7  Normalized  transition  kernel  versus  energy  transfer  for 
V  -  3.16  and  p  »  .5  with  C  -  10  and  30.  Similar  results 
are  obtained  for  all  p  and  v. 

8  Energy  transfer  moments  versus  incident  energy  for  v  -  3.16 

and  u  -  .5.  Criterion  ^  -  1/2  (3A2/3e)  is  satisfied 
for  C  >  2. 

9  Plot  of  mean  square  energy  transfer  per  collision  versus 

incident  energy  for  V  -  3.18,  p  »  .72  and  v  -  3.1*,  p  -  .46. 
Iteration  is  useful  for  A2MQ  >  'v-  1  and  the  diffusion  equa¬ 
tion  is  useful  for  <  *  1. 

10  Steady  state  reduced  adatom  energy  population,  (a)  diffusion 
approximation  for  v  «  3.16  and  p  -  .5  with  6-10  and  30; 

(b)  comparison  of  iteration  and  diffusion  approximations 
for  V  -  3.18  and  p  -  .72  at  6  -  55,6.  X  indicates  second 
Iteration,  x2W). 

11  Ratio  of  the  steady  state  to  the  equilibrium  rate  constant 
versus  frequency  parameter,  v  for  6-10  with  contours  of 

p  -  .2  and  1  shown.  The  steady  state  rate  constant  is  orders 
of  magnitude  below  the  equilibrium  rate  constant  at  "resonance" 
values  of  v. 

12  Resonance  values  \>r(p)  of  the  frequency  parameter  in  p, 
v  space.  The  adiabatic  limit  is  reached,  and  the  mean 
energy  transfer  goes  to  zero  on  the  contours  Indexed  by  n. 


Nondlmenslonal  steady  state  end  equilibrium  adsorption 

rate  constants  versus  veil  depth  for  v  ■  3.18,  y  -  .72 

and  v  ■  2.66,  y  -  .38. 
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H(0)  n  versus  the  frequency  parameter  v  for  the  Inertia 

parameter,  y  •  .2  and  1  with  the  temperature  parameter, 

q  ■  0,  1,  and  10.  At  resonance  and  at  n  •  1,  11(0)  »  1. 
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H(0)  n  versus  the  ratio  of  the  surface  temperature  to 
the  gas  temperature,  n  ■  T/T  for  v  •  3.18.  y  -  ,72,  and 
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